(N 

o 

(N 



(N 



Quantum State Transfer in a Two-dimensional Regular Spin Lattice of Triangular 

Shape 

Hiroshi Miki and Satoshi Tsujimotcu 

Department of Applied Mathematics and Physics, Graduate School of Informatics, 

Kyoto University, Sakyo-Ku, Kyoto 606 8501, Japan 

Luc Vine10 

Centre de recherches mathematiques, Universite de Montreal, 
P. O. Box 6128, Centre-ville Station, Montreal (Quebec), H3C 331, Canada 

Alexei Zhedano\[3 

Donetsk Institute for Physics and Technology, Donetsk 83 114, Ukraine 

(Dated: February 25, 2013) 

Quantum state transfer in a triangular domain of a two-dimensional, equally-spaced, spin lat- 
tice with non-homogeneous nearest-neighbor couplings is analyzed. An exact solution of the one- 
excitation dynamics is provided in terms of 2-variable Krawtchouk orthogonal polynomials that 
have been recently defined. The probability amplitude for an excitation to transit from one site to 
another is given. For some values of the parameters, perfect transfer is shown to take place from 
the apex of the lattice to the boundary hypotenuse. 
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The purpose of this letter is to present an exact so- 
lution of the 1-excitation dynamics of a regular two- 
dimensional spin-lattice of triangular shape with inhomo- 
geneous nearest-neighbor interactions. It also aims to ex- 
amine the quantum state transfer properties of this model 
in the plane. The analytic treatment is proving possible 
because of recent advances in the theory of multi- variable 
orthogonal polynomials yj, [2| . The transfer of quantum 
states between locations is fundamental in quantum in- 
formation processing. (See [3, |j] for a review.) Two 
desirable features are that the transfer be perfect, i.e. 
realized with probability 1, and that no external control 
be required. It is known that spin chains can model such 
one-dimensional perfect wires with the intrinsic dynam- 
ics generating on its own the transfer between the input 
and output sites with probability 1 after some time. 

The paradigm example 5] is that of the XX chain with 
Hamiltonian 



JV-l 



H = \ E ■ft+iW °f+i 



'i u i+ 



i) 



(i) 



1=0 



preserve the number of excitations and many aspects of 
quantum state transfer can be studied by restricting to 
the single-excitation eigenspace. In this case, the Hamil- 
tonian takes the form of a tri-diagonal (Jacobi) matrix 
and it can hence be diagonalized by orthogonal polyno- 
mials. 

The construction of XX spin chains with nearest- 
neighbor couplings that allow perfect state transfer 
(PST) is by now well understood 0,Q as an inverse spec- 
tral problem. For PST to occur, the 1-excitation spec- 
trum must be such that successive eigenvalues x s , s — 
0, !,■•• ,N satisfy 



Xs+l 



= ^M s 



(3) 



where T is the time of transfer and where M s are arbi- 
trary positive odd integers. Once a spectrum satisfying 
this condition is given, there is an algorithm [6| to ob- 
tain the (non-homogeneous) coupling strengths between 
sites that determine the XX Hamiltonian. The simplest 
spectrum that satisfies condition ([3]) is the linear one: 



with 



N 



s = 0, 1, 



,N. 



(4) 



J t = -y/l(N+l-l). 



(2) 



Here N + 1 is the number of sites and (jf,af (and a*) 
are the Pauli matrices acting at site I on a copy of C 2 . 
Such XX spin chains with nearest-neighbor couplings 
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The constructive approach of 6] directly shows that the 
corresponding couplings are those given in (J5|) which are 
recognized to be the recurrence coefficients of the sym- 
metric Krawtchouk polynomials J7J. We succeeded in 
finding the two-dimensional analogue of the Hamiltonian 
(JT]). We present it here and study its state transfer prop- 
erties. 

The letter is organized as follows. We first introduce 
the two-dimensional lattice Hamiltonian. It is defined on 
a triangular domain and contains 4 parameters in addi- 
tion to the integer N which specifies the number of sites 




the states 



FIG. 1. Two-dimensional regular lattice with site labeled by 
the couples (i,j), i,j G {0, ■ • ■ ,N} satisfying the triangular 
condition i + j < N. 



\i,j) = E i , j , i,j = 0,---,N 



(7) 



where Eij is the (N + 1) x (N + 1) dimensional matrix 
with a 1 at the (i,j) entry and zeros everywhere else. 
Note in relation with the action of the Pauli matrices, 
that the 1 and the in this register stand for the vectors 
( ) and (J respectively of the corresponding sites. We 
shall now look for the one-excitation eigenstates of H 
denoted by \s,t) with eigenvalues x s y. 



H\s,t) = x a ,t \s,t) . 



(8) 



(^). We then provide its one-excitation eigenstates in 
terms of two- variable Krawtchouk polynomials. This al- 
lows to examine transition amplitude for an excitation 
to move from one state to another. Particular atten- 
tion is paid to the case where the walk is initiated at 
the apex of the triangle. For some special values of the 
parameters, the probability to reach a site of the bound- 
ary hypotenuse follows a binomial distribution. Perfect 
transfer is thus seen to occur between the apex and the 
set of lattice points on that boundary. 

Consider the two-dimensional spin lattice depicted in 
fig. [TJ Its points are labelled by two indices (i, j), i,j € 
{0, • • • , N} such that i + j < N. The Hamiltonian will 
be taken of the form 



*= E 



J-i+l.j < x x 



KX+lj+<^+l,j) 



0<i+j<N 



with Iq j = Ji, 



Jj,j+1 , x x , v y \ 
2 V <7 i,i cr j,j+i + <J ij cr i,j+i) 



and I{ j = = Ji 



(5) 



constants 1+j and Ji j are couplings the sites (i 



,j i + j > N. The 
l,j) and 
(i,j) and the sites (i, j — 1) and (i,j) respectively. In the 
last term, Bij are the strengths of the magnetic fields 
at the sites (i,j). Such Hamiltonians are invariant under 
rotations about the z axis and thus preserve the total 
number of spins over the lattice that are up (or down): 



0<i+j<N 



(6) 



Introduce the expansion of \s, t) in terms of the \i,j) ba- 
sis: 



\s,t)= Y, Wu(a,t)\i,j). 

0<i+j<N 



(9) 



The eigenvalue equation ([8]) is seen to impose that the 
coefficients Wij(s,t) satisfy the 5-term recurrence rela- 
tion 

x s ,tWij= I i+ x t jWi + ij + Ji,j+iWij + i 

+B l . J W l . J + IijWi-i d + JijWij-i. (10) 

Let us now fix the coupling constants and the magnetic 
fields in ([5]) to be given by the following expressions which 
involve the four real parameters p\,p2,P3 and p<± and N: 

1 



Pi 



P?, 

1 



-\/Sp 1 p 3 {p2 +p 4 )i(N + 1 -i-j), 
V s P2Pi(Pi + Ps)j(N + 1 - i - j), 



B iJ = 



P2 +P4 

(N-i- j)S ( P2P4 (P1+P3) PlP3 (P2 + Pi ) 
PlPi~P2P3 I 



+3 



PiPi - P2P3 

P2 +P4 



P2 +Pi 

.PlPi - P2P3 
I 

Pi +P3 



Pi +P3 



(11) 



with S = p\ + P2 + P3 + Pa- With these coupling con- 
stants and magnetic fields, an exact solution is obtained 
in terms of 2- variable orthogonal polynomials. Precisely 
let 



W i>j (s,t) = W , (s,t) 



Kjj(s,t) 



with 



(P1P4-P2P3) 2(1+3) S-^ fN 

{PlPz{P2 +P4))' l {p2Pi(Pl +P3)) j \i,j 



(12) 



(13) 



We shall now restrict H to the subspace spanned by the 
states which contain only one excitation, that is only one 
spin up. A natural basis for that subspace is given by 



where ( . .) — .. ..,,. '. — ^r is the trinomial distribution 
function. Eq. (|10[) now amounts to the following relation 
for Kij{s,t) 



r- , , x ,1 PlP3(P2+P4)S P2P4JP1 +P3)S 



(Pi +P3)(PlP4-P2P3) 



(P2 + Pi){PlPi - P2P3) 



■(K itj+1 (s,t) - K iti (a,t)) }+i PlP4 ■ P2P3 (Jr i _ 1|J -(a,t)-Jr«, 3 -( 8 ,t))-j PlP4 ■ P2P3 (i^-i(M) - #<>,*))• 



Pi +P3 



P2 +P4 



I 

This is recognized as the 5-term recurrence satisfied by a linear expression in s and i. These 2-variable 
the 2-variable Krawtchouk orthogonal polynomials that Krawtchouk polynomials have the following explicit ex- 
have been introduced and characterized recently [l|, |2|. pression: 
The spectrum is obtained from this identification and 
thus found to be 



x s ,t = (pi+P2)s-(p 3 +P4:)t, 0<s + t<N, (15) 



K itj (s,t) - 



E 



(-j)k+l{- j)m+n(- S)k+m{-t)l+n i j k I 



0<k+l+m+n<N 



k\llm\nl(-N) 



k-\-l-\-m-\-n 



L 1 U 1 U, 2 U 2 



(16) 



with 



{Pi + P2){Pl + P3) {Pi + P2)(P2 + P4 

ui= 3 ,ua = 

(Pi +P3)(P3+P4) 
«1= 7, > V 2 = 



P2S 
(P2 +P4)(P3+P4) 



P35" 



P4S" 



(17) 



and with (a)^ = a(a + 1) • • • (a + fi — 1) the standard 
Pochhammer symbol. As eigenstates of a Hcrmitian 
Hamiltonian, the vectors \s, t) can be taken orthonormal: 

(a,t\u,v) = 5 SlU 6 t>v , (18) 

while the one-excitation vectors manifestly are 

(ii,ji\i2,J2) = 6i u i2^ji,h- ( 19 ) 

This requires the transformation matrix Wij{s,t) be- 
tween these basis to be orthogonal. Knowing the weight 
for which the polynomials Kij(s,t) are orthogonal, we 
find that this imposes that 



W ,o(a,t) 



N\ 1 

S,t)[{ Pl +P3)(P2+P4)} N 



(PlPi "P2P3) 2(JV S t KPlP2) S (P3P4Y 
(pi+P2) N - t (p3+P4) N - S 



(20) 



This completely determines the elements Wi,j(s,t). As 
a companion to ([9]). in view of the orthogonality of 
Wi ! j(a,t), we have the inverse relation 



l*.i)= E Wi,j(a,t)\a,t). 

0<s+t<N 



(21) 



We can now write down an explicit expression for the 
transition amplitude ftij),(k,l) (^) f° r an excitation at the 



site (i,j) to be found at the site (k, I) after some time T: 
/(«),(M)( T ) = (»,i|exp(-»Tir)|fc,0 (22) 

= E ^j(s,i)VK fc , ; (s,i)e lT[(p3+P4)t " (pi+P2)s] . 

0<s+t<iV 

(See 3 for analogous calculations in the case of quantum 
chains.) From now on we shall restrict the parameters 
and consider the case where 



Pl=P4, P2=P3- 



(23) 



With z = e -' T (pi+P2) ; the transition amplitude ||22J| sim- 
plifies under (|2"3")l to 

f(iJ),{k,l)( T ) = E Wij(*. ^W'fc.K*. *>*"*• ( 24 ) 

0<s+t<iV 

These quantities can be calculated with the help of sum- 
mation formulas 1]. In particular, if the quantum walk 
is initiated at the site (0, 0) we have 



f(p,o),(i,j)(T) - 



(2 Pl p 2 (z-iy + (p 1 +p 2 yz) 



2 7 \N-i-j 



(Pi - P2) l {z - l) i+j (p2Z + p 1 ) i (p 1 Z+P2) j 



Now take 



(P2- P i)-iz N ( Pl +p 2 y N 



T 



Pi +P2 

(i.e. z = —1), the amplitude (|2"5")l becomes 



(25) 



(26) 



f (0,0), (i,j) 

\Pl ~ P2 
(Pl+P2) 2N 



TT 



Pi +P2 

(i+3) 



'(2P1P2, 



i+j 



N 
i,3 



&P1P2 ~ (Pi + P2 



2 \N-i-j 



(27) 



Imposing further that 

(pi +P2) 2 = 8P1P2, 
one arrives the following simple expression 

f(0,0),(i,j) 



Pi +P2 




i+j,N- 



(28) 



(29) 



The output excitation will hence distribute binomially on 
the site of the boundary hypotenuse. It is moreover easy 
to see therefore that 



N 

E 

k=0 



f(0,0),(k,N-k) 



Pi +P2 



(30) 



We thus observe that with these values of the parame- 
ters, the Hamiltonian will evolve the state 1 0, 0) in time 

T = T„ to any one of the states \i, N — i) with proba- 



P1+P2 



bility one. This may be looked at as some kind of perfect 
transfer. As a consequence we have that 



/(0,0),(t,; 



,Pl +P2, 

This is observed to generalize to 



= 0, i + j <N. 



(31) 



f(i,j),(k,i) 



Pi +P2 



= 0, i + j + k + I < N. (32) 



It is interesting to see what is the form of the Hamil- 
tonian when the choices (|2"3"f and (|2"5|) are made for the 
parameters. One finds that the coupling constants Iij 
and Jij and the strengths Bij of the magnetic field are 
given by the following simple expression: 

1 



*i,j 


Pi +P2 


Pi +P2 
B i,j 



= -^Vj(n + i-i-j), 



Pi -Pi 



Pi +P2 Pi+ P2 



U-i)- 



(33) 



Here it should be noted that (|28| implies that 



Pi - P2 
Pi +P2 



± 



V2' 



(34) 



We see that the magnetic fields are antisymmetric with 
respect to the lattice diagonal and we observe the symme- 
try Iij —Ji t j under % •<-> j. When this specialization of 
the parameters is examined at the level of the 2-variable 
Krawtchouk polynomials, symmetric functions in s and 
t result. 

Let us summarize to conclude. We have presented 
a spin lattice which is the two-dimensional analogue of 
the simplest non-homogeneous quantum spin chain that 
exhibit perfect state transfer. In parallel to the one- 
dimensional case, the one-excitation Hamiltonian is diag- 
onalized by two-variable Krawtchouk polynomials. This 
exact solution allows to show that there is some kind 
of (generalized) perfect state transfer between one input 
site and a set of output sites. One can imagine various 
ways in which this could prove useful. Very few spin 
lattices are known to be analytically solvable [9]. Al- 
though we only have provided an exact determination of 
the one-excitation dynamics, we trust that this addition 
will allow to further explore phenomena inherent to two 
(higher) dimensional spin networks. 
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